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Abstract 

Central Pattern Generators (CPGs) in invertebrates are comprised of networks of neurons in 
which every neuron has reciprocal connections to other members of the CFG. This is a "closed" 
network topology. An "open" topology, where one or more neurons receives input but does not 
send output to other member neurons, is not found in these CPGs. In this paper we investigate 
a possible reason for this topological structure using the ability to perform a biological functional 
task as a measure of the efficacy of the network. When the CFG is composed of model neu- 
rons which exhibit regular membrane voltage oscillations, open topologies are essentially as able 
to maximize this functionality as closed topologies. When we replace these models by neurons 
which exhibit chaotic membrane voltage oscillations, the functional criterion selects closed topolo- 
gies when the demands of the task are increased, and these are the topologies observed in known 
CPG networks. As isolated neurons from invertebrate CPGs are known in some cases to undergo 
chaotic oscillations ||Hayashi fc Ishizuka, 1992| , [Abarbanel et al 1996|| this provides a biological basis 



for understanding the class of closed network topologies we observe. 



1 Introduction 

Central Pattern Generator (CPG) neural networks in invertebrates perform a wide variety of func- 
tional roles, each of them requiring rhythmic output from the CPG to the muscles controlling the 
function [[Marder fc Calabrese 1996| . In our study of the pyloric CPG of the California spiny lobster 
Panrulis interruptus we have addressed the question why the network is closed in the sense that 
all neurons in the CPG both send signals to other neurons in the circuit and receive input from 
other neurons in the circuit. None of the neurons only receive from or send signals to the rest of 
the network. The pyloric CPG is illustrated in the right half of Figure lA. 

In this paper we address the appearance of this closed CPG network and seek a functional 
biological reason why evolution may have been led to select this configuration. We are aware of 
investigations [[Getting 1989|| which attempt to answer this question in terms of enhancing robust- 



ness or self-organization. In the present work we test an alternate hypothesis: the closed network 
topology composed of irregular neurons provides greater efficiency in the task which the network is 
required to perform. We present here evidence for this hypothesis in terms of calculations of the 
ability of a model crustacean pyloric CPG to perform its task of transporting shredded food from 
the stomach to the digestive system. We build the CPG first from conductance based Hodgkin- 
Huxley (HH) [Podgkin fc Huxley, 1952] , Puerta 1996|| neurons which exhibit periodic oscillations 



and subsequently from model neurons of Hindmarsh-Rose (HR) ||Hindmarsh fc Rose, 1984| type 



which have chaotic membrane voltage oscillations. Using the transport criterion made explicit be- 
low, we find the periodic HH neurons to allow open topologies in CPG operation while chaotic 
component neurons allow only closed topologies on our functional criterion. 

In an earlier paper fHuerta et al, in press|| we described in some detail the use of regular HH 



neurons with two spatial compartments as CPG network components. Many of the allowed network 
configurations, were open and quite robust in the presence of noise as the HH models oscillated 
in a strongly dissipative limit cycle and did not sit in parameter space near a bifurcation point. 
In real biological networks the component neurons show substantial irregularity in their bursting 



cycles |[Elson et al 1999| , [Elson et al 1998|| , and this has led us to the work reported here. 

The chaotic model neuron we use in our investigations is a three dimensional version of an HR 
neuron. It exhibits spiking-bursting behavior as observed in the laboratory, and this alone would 
not distinguish it from similar behavior of HH type neuron models. Some evidence, peripheral to 
this study, for the realistic structure of the HR neurons comes from its observed ability to act as 
a substitute for biological neurons in the lobster pyloric CPG when it is realized in simple analog 
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Figure 1: Examples of closed topolog ies in invertebrate CPGs: (A) the gastric and pyloric C PGs in Crustacea (mo d- 
ified from [ ^elverston fc Moulins 1987 |), (B) The f eeding CPG in Planorbis (modified from [ Arshavsky et al 1985| ), 
(C) The swimming CP G in Tritonia (modified from [ Getting 1989t ), (D) The swimming CPG in Clione (modified from 



I Arshavsky et al 1991 . Dots represent inhibitory synapses while arrows represent excitatory synapses. Electrotonic 
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gap junctions are represented by resistors 



circuitry and coupled into that biological network ||Szucs et al 1999[ . The main feature of the chaotic 
neuron is that it is built around a homoclinic loop, which carries the fast oscillations. Motion on the 
homoclinic loop is intrinsically unstable which means that any perturbation close to the homoclinic 



structure will produce striking differences in its membrane voltage time course: the neuron will 
either fire another spike or hyperpolarize the membrane [ Bazhenov et al 1998| , Abarbanel et al 1996 , 



Wang 1993 |. The dynamical pictures of the chaotic and the regular models are completely different. 

The functional criterion we investigate is based on the role of the pyloric CPG in the California 
spiny lobster. The pyloric CPG has as its task the transport of shredded food from the stomach to 
the digestive system. Our functional criterion for the selection of network topologies is the ability to 
perform this transport as the severity of the task is increased. In an intuitive way one might expect 
to maximize this transport through regular CPG oscillations, and this is consistent with a fully 
coupled or closed topology of interconnections which regularize the chaotic behavior of individual 



components of the CPG assembly pazhenov et al 1998|| . Our overall picture then suggests a cor 



respondence between the functional task of a CPG, the observations of regular patterned behavior 
in CPG assemblies of individually chaotic components, and the CPG network topology. 

2 The chaotic model neuron 

In our simulations we used a modified, three dimensional version of the HR model neurons. The 
model is comprised three dynamical variables comprising a fast subset, x(t) and y(t), and a slower 
z{t). x{t) represents the cell's membrane potential. These dynamical variables satisfy 



dx{t) 
dt 

dy{t) 
dt 
1 dz{t) 



%(t) + 1.5x2(t) - 0.25x3(t) - 2z{t) + 2e + Jsyn (1) 

0.5-0.625x^(t)-y(t), (2) 

^^ -z(t)+2[x(t) + 3.2] (3) 

where e represents an injected DC current, and yU is the parameter that controls the time constant 
of the slow variable. The parameters were chosen to place the isolated neurons in the chaotic 
spiking-bursting regime: e = 3.281, ^ = 0.0021. 

Isyn represents the postsynaptic current evoked after the stimulation of a chemical graded 
synapse. In this paper we consider only inhibitory synapses, which is the main kind of inter- 
connection present in the pyloric CPG of the lobster. The synaptic current has been simulated with 
the traditional description used in the dynamical clamp technique [pharp et al 1993| with minor 
modifications: 

%n = -^'^(a;pre)^(a;post) (4) 



where g is the maximal synaptic conductance, and XpQg^^ is the membrane potential of the postsy- 
naptic neuron. r[xpre) is the synaptic activation variable determined from the presynaptic activity 
by: 



7 , [' oo \-^pre ) '\l'r 

roo{Xpre) = [l+taiih{{xpre + A)/B)]/2. 



(5) 
(6) 



In our work A = 1,2 and B = 0.9. r^. is the characteristic time constant of the synapse (x^ ~ 100). 
'^{xpost) is a nonlinear function of the membrane potential of the postsynaptic neuron: 



'^yxpgsfj 



[1 + tanh((a;post + o,)/b)). 



(7) 
(8) 



The parameters for "^{xpost) (« = 2.807514 and b = 0.4 in this paper) were chosen so that this 
function remained linear for the slow oscillations, the subthreshold regions for the fast spikes: 



d^{x) 



dx 

i^{x) 



a,b 



a,b 



0. 



(9) 

(10) 



Each pair of model neurons were connected in mutual inhibitory configurations consisting of 
basic two-cell units such as the one shown in Figure 0. 




Figure 2: Basic configuration of synaptic couplings of the model for mutual inhibition between two neurons. The 
graded synapses are described by equations (B]|) in the text. 



3 The models for the CPG and the mechanical device 



We have previously built a CPG that controls a pyloric chamber model jHuerta et al, in press] 
The CPG is composed of three neurons with mutual inhibitory coupling as shown in Fig. 2 



6 



The configuration we selected is shown in Figure 3. The six independent maximal conductances 
9ij'^ i 7^ J5 hi = 1)2,3 are to be selected as described below. The CPG sends electrical activity to 
the 'muscles' which control the dilation and contraction of the 'pyloric chamber' represented here 
by a simulated mechanical 'plant'. Three neurons are the minimal number of neurons that produces 
a maximization of the average flow of food pumped out of the plant. 
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Figure 3: Illustration of the mechanical plant and the CPG. In the upper part of this figure we show 
the three neuron CPG and the six possible connections Qij among them using mutual inhibition. In 
the lower half of the figure we show the pumping 'plant' composed of three pipes which scan vary 
their radius. Eac h CPG neuron only aff ects one segment of the pipe. A detailed description of the 
pump is given in IPuerta et al, in press| . 



The mechanical plant is inspired by the lobster pyloric chamber ||Selverston fc Moulins 1987[| . In 
Fig. 3 we have sketched the model with three joined pipes that can expand or contract radially. The 
central pipe only dilates and the end pipes only contract; this mimics the behavior of the pyloric 
chamber of the lobster. The proper dynamical combination of muscle activity influencing the walls 
of the pipes will lead to movement of shredded food towards the right end of the pyloric chamber. 

In our selection procedure for combinations of the g^j, we set the average transport through the 
pump in time T to a set of increasingly larger values and determine which configurations are able 
to achieve each level of transport. The average throughput over time T is given by 

fto+T 

Hg^,) = ^ / As{t)v,{t)dt, (11) 



to 



where p is the density of the material being pumped, A^{t) is the cross sectional area of the 
rightmost pipe section, VQ{t) is the mean velocity through the rightmost cross section of the plant, 
and to is some starting time. The differential equations determining A:i{t) and VQ{t) were derived 
Huerta et al, in pressj] . They come from the Navier-Stokes equations and from mass and energy 



m 



conservation. We treat the shredded food as homogeneous, incompressible and isothermal, and it 
moves with low Reynolds number so the flow is always laminar and all radial velocities in the pipe 
sections are small. We also assume that no food leaks out during the pumping and that there were 
no head-losses in the joints of the pipes. 

In our calculations the mean initial velocities in the sections are zero. We choose a set of gij as 
indicated in a moment, and ask whether the specified configuration can achieve a chosen level of 
$((7jj). The system of differential equations was integrated using a Runge-Kutta 6(5) scheme with 
variable time step and with an absolute error of 10~^^ and a relative error of 10~^. 

We choose the gij to be 0, 50nS, or 200nS, and evaluate ^{gij) for each of these 3^ choices. The 
system was run for a simulated time of 60 sec; a transient of to = 10 sec was eliminated and ^{gij) 
was found by averaging over the last T = 50 sec. 

In the following section we will compare the performance of both regular and chaotic neurons 
in controlling the pump for the different connection architectures. The detailed description of the 
regular HH model can be found in [[Huerta et al, in press||. The chaotic neurons were implemented 



as described in the previous section. Both types of neurons are able to generate spiking-bursting 
behavior. 

4 Analysis of solutions 

< / CPG j 

Neuron ^^^^^^y/^ 

Figure 4: Definition of an open network topology: there must be one or more neuron without 
feedback from the remaining neurons. 

It is useful to distinguish among three different types of topologies: 




• Open topology: in this case there is one or more neuron which receives no feedback from 
the other network elements. This is shown in Figure 4. 

• Semi-open topology: in this case there is one or more neuron which receives input from 
some neurons in the network, but no connection is made back to the network. 

• Closed topology: in this case there are connections from each network member to and from 
the remainder of the network. 

In order to quantify the difference between open topologies and the others we introduce the 
following definitions. We define the set of solutions with a resulting flow lying in ^lower < *& < ^upper 
as 

G{^lower, ^upper) = {{9l2, 913, 921, 923, 931, 932)', ^lower < $ < ^upper} ■ (12) 

and the subset of open solutions as 

Oi^lower, Supper) = {{9l2, 9l3, 921,923, 931,932)] ^i/9ji = 0, Vj y^ V, ^lower < $ < ^ upper} ■ (13) 

We report the ratio between the number of open solutions | 0{^iower, ^upper) \ and the total number 

of solutions I G{^lower, Supper) \ 



/I I \ 1 I '^y^lower, ^upper) \ / | ^y^lower, ^upper) \ '^ \ ^\^lower, ^upper) \t^ U /-. .\ 

Vy^lower, ^upper) ~ \ f^ .r i /^/^ ^ \ I n \ ) 

U II I <~J [(^ lower, ^ upper) \— U 

as a quantitative measure of the topology of allowed solutions at a given value of $. | . . . | denotes 
the cardinality of a set. Since | 0{^iower, ^upper) \<\ G{^iower, ^upper) |,0 < r/ < 1. 
We investigated two questions in the calculations reported here: 

• Does the use of chaotic component CPG neurons filter out open solutions ? 

• Can noise eliminate open solutions when the component CPG neurons are regular oscillators? 

In examining the first question we used three different values of the Qij as indicated above for 
both the regular neuron model and the chaotic neuron model. The number of configurations was 
3^. The number of simulation trials for a given configuration pattern was 15 for different random 
initial conditions. In Figure || we can see a plot of '7('^iower' "^upper) as defined in Equation [14. 



We can see that as the fiow $ increases there are no open topologies of chaotic component neurons 
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which can achieve the flow level. $ must be substantially reduced to find allowed open topology 
synaptic connection configurations. 

If we build the CPG with regular neurons we see that there are many open topology solutions 
for high values of $. Moreover, the ^($) is strongly decreasing for the CPG with chaotic neurons. 
However, the intrinsically regular CPG does not lead to a strong dependence of r] on the level of 
the fiow. 

This is one of the essential points of our work. Since the neurons in the pyloric CPG appear to 
oscillate chaotically, a closed topology of synaptic connections is required to achieve high transport 
rates for food. Regular neurons can achieve high transport rates in open and in closed configurations, 
but the pyloric CPG is observed to be closed. Our calculations provide a suggestive connection 
between these observations. CPGs made of chaotic neurons could have evolved open topologies, 
but these are not as effective in transporting food as closed topologies and presumably were not 
selected in evolutionary processes. This means that by means of chaotic neurons a selective choice 
of available solutions is obtained. Intrinsic chaotic dynamics in the isolated neurons provide a 
straightforward explanation of the existence of non-open connection topologies in invertebrate's 
CPGs. 

To address our second question iid N(0,1) white noise e{t) is introduced in the external injected 
currents of the Hodgkin-Huxley model Puerta et al, in press|| so I(t) = lo + creit) where a is the 
amplitude of the noise. We introduce this type of additive noise for simplicity, because it can be 
more easily integrated [ [Mannella fc Palleschi[| since there is no dependence on the other variables 



of the ordinary differential equations. In Figure 5 we can see the results of our calculations for 8 
different values of a. For each a value 20 trials of different initial conditions for a given configuration 
of Qij were carried out and the average of those values was taken as the output of the calculation. 
The quantity we used to estimate changes as a function of a is 

where ^max = maxg{$} and ^th is the threshold value that determines the existence of good 
solutions. If a reduction in the number of open configurations is observed the value 77(0") must 
decrease. If our selected configurations are eliminated by noise, then open-topology configurations 
should disappear and f]{cr) must tend to 0. Otherwise for any value of noise then r]{a) is always 
greater than 0. In fig 5 we can clearly see that for any value of the noise no reduction of the 
open-topology configurations is observed. Moreover, the quantity rj approaches 1 for high values 
of the noise which answers our question about noise in the positive: the open configurations of 
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Figure 5: A. ''^('^lower' "^upper) for the regular and chaotic component CPGs. The data are pre- 
sented by placing the average transport in bins defined by "^lower ~ ^ ^^^ ^upper = 1.01$. As 
we increase the transport level $, the number of open solutions which can achieve this level goes 
rapidly to zero. B. Number of total allowed configurations for both the regular and chaotic CPG 
as a function of the level of average transport $. 
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regular component CPGs are robust against the presence of synaptic noise. The appearance of 
closed configurations must have a different explanation, and that we have provided above. 
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Figure 6: 77(a) for regular Hodgkin-Huxley neurons as CPG components, model. The open-topology 
configurations are not eliminated because of the noise. In these calculations $ = 0.06 (gr/sec). 



5 Conclusion 



In this paper we have concentrated on the lobster pyloric CPG as we have experimental evidence of 
the chaotic oscillations of its component neurons when they are isolated from the intact network. 
There are many other invertebrate CPGs as shown in Figure 1, and a common thread among these 
CPGs, whether they are involved in shredding food, as in the gastric CPG of Crustacea, feeding, 
as in the Planorbis, or swimming as in Tritonia is that the topology of the CPG network is closed. 
It would be very interesting to follow the path we have set out here for these and other CPGs to 
establish whether the connection between network topology , efficiency in biological function and 
chaotic neuronal elements holds there as well. 



Since the early papers that showed evidence of chaos in CPG neurons [[Hayashi fc Ishizuka, 1994 
[Mpitsos et al 1988| little attention has been given to the role of chaotic neurons in CPGs. It 
appears that the main reason is that intact CPG networks work in a regular fashion. Indeed 
their task is to provide a regular rhythm to the animal so that the appropriate response of the 
muscles that operate in a physical device, such as the pyloric chamber, performs a specific function. 



12 



Therefore, the interest of modelers was focused mainly on studying the dynamics of periodic model 
neurons [[Skinner et al 1993| , [Kepler et al 1990| , [lerman et al 1998| , [Baer et al 1995[ , [Roberts 1998[ , 



Huerta 1996[| . In very few cases modelers have tried to understand the potential advantages of 
using chaotic neurons [ [b'reeman 1996| , [Kabinovich fc Abarbanel 1998[ , [Zak 199l| , but typically in 
more complex systems, not in CPGs. By means of this work we hope to draw the reader's attention 
to an interesting set of roles played by chaotic neurons in CPGs. 

The connectivity of neurons in CPGs is probably determined by many developmental and evo- 
lutionary factors, and in this paper we have considered an aspect of the required function of a CPG 
as providing a significant driving force in this topological decision. Starting from the observation 
that component neurons from the lobster pyloric CPG have chaotic membrane potential oscillations 



when observed in isolation [[Abarbanel et al 1996[| , we have investigated the influence that chaotic 
versus regular neurons might have on the structure of the neural interconnections in the CPG in 
order to achieve its functional goal. If our conclusions are to extend to other CPGs, clearly one 
must systematically, as we have for the pyloric CPG of lobster, establish the modes of oscillation 
for the CPG members in isolation. We conjecture that the example of lobster pyloric CPG is not 
special, but we are seeing in this example a pattern which will repeat in other CPG systems. 
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